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Abstract 

We review some old and prove some new results on the survival probability of a random walk 
among a Poisson system of moving traps on 1 d , which can also be interpreted as the solution of a 
parabolic Anderson model with a random timc-dcpcndcnt potential. We show that the annealed 



survival probability decays asymptotically as e Aiv/ * for d = 1, as e A2 '/ log * for d = 2, and as 



-Xdt 



for d > 3, where Ai and A2 can be identified explicitly. In addition, we show that the quenched 
(T) ■ survival probability decays asymptotically as e~ Xdt , with > for all d > 1. A key ingredient in 

bounding the annealed survival probability is what is known in the physics literature as the Pascal 
principle, which asserts that the annealed survival probability is maximized if the random walk 
stays at a fixed position. A corollary of independent interest is that the expected cardinality of the 
range of a continuous time symmetric random walk increases under perturbation by a deterministic 
<^ path. 

AMS 2010 subject classification: 60K37, 60K35, 82C22. 

Keywords: parabolic Anderson model, Pascal principle, random walk in random potential, trapping 
dynamics. 
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1.1 Model and results 

Let X := (X(t))t>o be a simple symmetric random walk on % d with jump rate k > 0, and let 
(yf)i<j<N y ,y&z d be a collection of independent simple symmetric random walks on Z d with jump rate 
p > 0, where N y is the number of walks that start at each y £ Z d at time 0, (N y ) y& %d are i.i.d. Poisson 
distributed with mean v > 0, and Yj := (yJ(t))t>o denotes the j-th walk starting at y at time 0. Let 
?-h I us denote the number of walks Y at position x G Z d at time t > by 

C(t,x):= $»C*7(t)). (1) 



It is easy to see that for each i > 0, (£(t, x)) x ^d are i.i.d. Poisson distributed with mean v, so that 
(^(t, -))t>0 is a stationary process, and furthermore it is reversible in the sense that -))o<t<T is 
equally distributed with (£(T — i, -))o<t<T- We will interpret the collection of walks Y as traps, and 
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at each time t, the walk X is killed with rate "f£(t, X(t)) for some parameter 7 > 0. Conditional on 
the realization of the field of traps £, the probability that the walk X survives by time t is given by 

rt 



exp{- 7 jf £( S ,X(s))da}], (2) 



where Kq denotes expectation with respect to X with X(0) = 0. We call this the quenched survival 
probability, which depends on the random medium £. When we furthermore average over £, which we 
denote by E*, we obtain the annealed survival probability 

Et[Z^} = E«Eff [exp { - 7 X{s)) ds}] . (3) 

We will study the long time behavior of the annealed and quenched survival probabilities, and in 
particular, identify their rate of decay and their dependence on the spatial dimension d and the 
parameters k, p, v and 7. 

Here are our main results on the decay rate of the annealed and quenched survival probabilities. 

Theorem 1.1 [Annealed survival probability] Assume that 7 6 (0,oo], k > 0, p > and v > 0, 

then 

f8pt 



exp 



{-^^(1 + 0(1))}, d=l, 



exp {-1/^^(1 + 0(1))}, d = 2, (4) 

exp{ - \d,<y,K,p,vtO- + d>3, 
where \d,y k,p,v depends on d, 7, k, p, v , and is called the annealed Lyapunov exponent. Furthermore, 
i,k,,p,v — ^d,-y,o,p,u — + ~ 7^ — ); where Gd{0) '■— Jq j>t(0)dt is £/ie Green function of a simple 

symmetric random walk on Z, d with jump rate 1 and transition kernel pt(-)- 

Note that in dimensions 1 and 2, the annealed survival probability decays sub-exponentially, and the 
pre-factor in front of the decay rate is surprisingly independent of 7 £ (0, 00] and k > 0. The key 
ingredient in the proof is what is known in the physics literature as the Pascal principle, which asserts 
that in Q, if we condition on the random walk trajectory X, then the annealed survival probability 
is maximized when X = 0. The discrete time version of the Pascal principle was proved by Moreau, 
Oshanin, Benichou and Coppey in |19|. I2Q| . We will include the proof for the reader's convenience. As 
a corollary of the Pascal principle, we will show in Corollary 12.11 that the expected cardinality of the 
range of a continuous time symmetric random walk increases under perturbation by a deterministic 
path. 

In contrast to the annealed case, the quenched survival probability always decays exponentially. 

Theorem 1.2 [Quenched survival probability] Assume that d > 1, 7>0, k > 0, p > and 
v > 0. Then there exists deterministic \d,i,n,p,v depending on d, 7, k, p, v, called the quenched Lyapunov 
exponent, such that F*-a.s., 

= exp{ - Xd,j,K,p,ut(l + o(l))} as t -> 00. (5) 

Furthermore, < \dj,K,p,is ^ l u + K f or all d > 1, 7 > 0, k > 0, p > and v > 0. 

Remark. When 7 < 0, Zj^ can be interpreted as the expected number of branching random walks 
in the catalytic medium £. See Section [1.31 for more discussion on this model. As will be outlined at 
the end of Section [4.11 ([5]) also holds in this case, and lies in the interval [—7^ — k, 00). 

In Proposition 13.21 below, we will also give an upper bound of the same order as in Theorem 11.11 
for the survival probability E^[Z^], where (£(0, x)) x& z d ^ s deterministic and satisfies some constraints. 
These constraints hold asymptotically a.s. for i.i.d. Poisson distributed (£(0, x)) xgZ d. Therefore we 
call this a semi-annealed bound, which we will use in Section [3] to obtain sub-exponential bounds on 
the quenched survival probability in dimensions 1 and 2. 
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1.2 Relation to the parabolic Anderson model 

The annealed and quenched survival probabilities and E^ [Z^t ] are closely related to the solution 
of the parabolic Anderson model (PAM), namely, the solution of the following parabolic equation with 
random potential £: 



d 

— u(t, x) = nAu(t, x) — 7£(i, x) u(t, x) 
u(0, x) = 1, 



xez d , t>o, 



(6) 



where 7, k and £ are as before, and A/(x) = ^ Yl\\y-x\\=i(f (v) ~ f( x )) * s the discrete Laplacian on 
Z d , which is also the generator of a simple symmetric random walk on Z rf with jump rate 1. 
By the Feynman-Kac formula, the solution u admits the representation 



x 



u(t,0) = Eq 



expj-7^ £(t-s,X(s))ds 



(7) 



which differs from Z^ in @ by a time reversal in £. When we average u(t, 0) over the random field 
£, by the reversibility of (£(t, -))o< s <t, we have 



E^[u(t, 



exp j-7^ £(t - s, X(s)) dsj = E^E« exp j- 7 ^ £(s, X(s)) ds 



E^]. (8) 



Therefore Theorem 1 1 . 1 1 also applies to the annealed solution E^[u(t, 0)]. Despite the difference between 
Z^£ and u(t, 0) due to time reversal, Theorem 1 1 . 2 1 also holds with u(t, 0) in place of Z^. 

Theorem 1.3 [Quenched solution of PAM] Let d > 1, 7 > 0, k > 0, p > 0, 1/ > and 

7,k,p,^ > oe i/ie same as in Theorem \LJ^ Then E^-a.s., 



t(t,0) = exp { - Xd niK)P)V t(l + o(l))} as t -»■ 00. 



(9) 



Remark. By Theorem 11.21 and the remark following it, for any 7 £ K, i" 1 logn(t, 0) converges in 
probability to —\d,^,K,p,v because ii(i,0) is equally distributed with Z t 7 ^. However, we were only able 
to strengthen this to almost sure convergence for the 7 > case, but not for 7 < 0. For a broader 
investigation of the case 7 < 0, see Gartner, den Hollander, and Maillard [T3], which is also contained 
in the present volume. 



1.3 Review of related results 

The study of trapping problems has a long history in the mathematics and physics literature. We 
review some models and results that are most relevant to our problem. 



1.3.1 Immobile Traps 

Extensive studies have been carried out for the case of immobile traps, i.e., p = and •) = £(0, •) 
for all t > 0. A continuum version is Brownian motion among Poissonian obstacles, where a ball 
of size 1 is placed and centered at each point of a mean density 1 homogeneous Poisson point pro- 
cess in Mr, acting as traps or obstacles, and an independent Brownian motion starts at the origin 
and is killed at rate 7 times the number of obstacles it is contained in. Using a large deviation 
principle for the Brownian motion occupation time measure, Donsker and Varadhan [7] showed that 

d 

the annealed survival probability decays asymptotically as exp{— Cd tl t d + 2 (1 + o(l))}. Using spec- 
tral techniques, Sznitman |24j later developed a coarse graining method, known as the method of 
enlargement of obstacles, to show that the quenched survival probability decays asymptotically as 
exp{— Cd,y ( log ^2/d (1 + Similar results have also been obtained for random walks among immo- 

bile Bernoulli traps (i.e. £(0, x) S {0,1}), see e.g. [HI SI [XI [2] - Traps with a more general form of the 
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trapping potential £ have also been studied in the context of the parabolic Anderson model (see e.g. 
Biskup and Konig [3]), where alternative techniques to the method of enlargement of obstacles were 
developed and the order of sub-exponential decay of the survival probabilities may vary depending on 
the distribution of £. Compared to our results in Theorems 11.11 and 11.21 we note that when the traps 
are moving, both the annealed and quenched survival probabilities decay faster than when the traps 
are immobile. The heuristic reason is that, the walk survives by finding large space-time regions void 
of traps, which are easily destroyed if the traps are moving. Another example is a Brownian motion 
among Poissonian obstacles where the obstacles move with a deterministic drift. It has been shown 
that the annealed and quenched survival probabilities decay exponentially if the drift is sufficiently 
large, see e.g. [Ml Thms. 5.4.7 and 5.4.9]. 

1.3.2 Mobile Traps 

The model we consider here has in fact been studied earlier by Redig in [22] . where he considered a 
trapping potential £ generated by a reversible Markov process, such as a Poisson system of random 
walks, or the symmetric exclusion process in equilibrium. Using spectral techniques applied to the 
process of moving traps viewed from the random walk, he established an exponentially decaying upper 
bound for the annealed survival probability when the empirical distribution of the trapping potential, 
7 Jo £( s ' 0) d s > satisfies a large deviation principle with scale t. This applies for instance to £ generated 
from either a Poisson system of independent random walks or the symmetric exclusion process in 
equilibrium, in dimensions d > 3. 

1.3.3 Annihilating Two- type Random Walks 

In [5], Bramson and Lebowitz studied a model from chemical physics, where there are two types of 
particles, As and Bs, both starting initially with an i.i.d. Poisson distribution on Z d with density 
Pa(0) resp. Pb(0). All particles perform independent simple symmetric random walk with jump rate 
1, particles of the same type do not interact, and when two particles of opposite types meet, they 
annihilate each other. This system models a chemical reaction A + B — > inert. It was shown in [S] 
that when pa(0) = Pb(0) > 0, then pA(t) and (the densities of the A and B particles at time 

t) decay with the order t~ d ^ in dimensions 1 < d < 4, and decay with the order t~ 1 in d > 4. When 
Pa(0) > Pb(0) > 0, it was shown that pA(t) — > /?a(0) — Pb(0) as t — > oo, and — log psit) increases 
with the order \ft in d = 1, t/logt in d = 2, and t in d > 3, which is the same as in Theorem 11.11 
Heuristically, as ps{t) — > and pA(t) — > Pa(0) — pb(0) > 0, we can effectively model the B particles as 
uncorrelated single random walks among a Poisson field of moving traps with density pa(0) — (0). In 
light of Theorem ll.il it is natural to conjecture that psit) decays exactly as prescribed in Theorem ll.il 
with v = pa(0) — Pb(0) and 7 = 00, whence we obtain not only the logarithmic order of decay as in 
[5], but also the constant pre-factor. However we will not address this issue here. 

1.3.4 Random Walk Among Moving Catalysts 

Instead of considering £ as a field of moving traps, we may consider it as a field of moving catalysts 
for a system of branching random walks which we call reactants. At time 0, a single reactant starts 
at the origin which undergoes branching. Independently, each reactant performs simple symmetric 
random walk on 7L d with jump rate k, and undergoes binary branching with rate |7|£(t, x) when the 
reactant is at position x at time t. This model was studied by Kesten and Sidoravicius in [15], and 
in the setting of the parabolic Anderson model, studied by Gartner and den Hollander in j!2j . For 
the catalytic model, 7 is negative in ([2]), ([3]), ([7]) and (jHJ), and and E^[Z t 7 ^] now represent the 
quenched, resp. annealed, expected number of reactants at time t. It was shown in [T2] that ^\Z]A 
grows double exponentially fast (i.e., t _1 log log [Z^] tends to a positive limit as t — > 00) for all 
7 < in dimensions d = 1 and 2. In d > 3, there exists a critical 7 Cj ^ < such that E^[Z^] grows 
double exponentially for 7 < 7^, and grows exponentially (i.e., t~ l log E^[Z t 7 £ ] tends to a positive 
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limit as t — > oo) for all 7 E (7c,<2i 0). In the quenched case, however, it was shown in [15] that only 
exhibits exponential growth (with log Zj^ shown to be of order t) regardless of the dimension d > 1 
and the strength of interaction 7 < 0. Such dimension dependence bears similarities with our results 
for the trap model in Theorems 11.11 and 11.21 

1.3.5 Directed Polymer in a Random Medium 

We used Z^ to denote the survival probability, because Z^ and E^[Zj^] are in fact the quenched 
resp. annealed partition functions of a directed polymer model in a random time-dependent potential 
£ at inverse temperature 7. The directed polymer is modeled by (X(s))o< s <t. In the polymer mea- 
sure, a trajectory (X(s))o< s <t is re- weighted by the survival probability of a random walk following 
that trajectory in the environment £. Namely, we define a change of measure on (X(s))o< s <t with 
density e" 7 ^ t(s,x( s ))ds in the quenc h e d model, and with density E^e" 7 ^ £( s ' X ( s )) ds ]/E^[^y in 
the annealed model. Qualitatively, the polymer measure favors trajectories which seek out space-time 
regions void of traps. However, a more quantitative geometric characterization as was carried out for 
the case of immobile traps (see e.g. |24j ) is still lacking. 

For readers interested in more background on the problem of a Brownian motion (or random walk) 
in time-independent potential, we refer to the book by Sznitman [23] on Brownian motion among 
Poissonian obstacles, and the survey by Gartner and Konig [11] on the parabolic Anderson model. For 
readers interested in more recent studies of a random walk in time-dependent catalytic environments, 
we refer to the survey by Gartner, den Hollander and Maillard |13j . For readers interested in more 
recent studies of the trapping problem in the physics literature, we refer to the papers of Moreau, 
Oshanin, Benichou and Coppey \19\ [20] and the references therein. 

After the completion of this paper, we learnt that the continuum analogue of our model, i.e., the 
study of the survival probability of a Brownian motion among a Poisson field of moving obstacles, 
have recently been carried out by Peres, Sinclair, Sousi, and Stauffer in [2"T] . See Theorems 1.1 and 
3.5 therein. 

1.4 Outline 

The rest of this paper is organized as follows. Section [2] is devoted to the proof of Theorem II .11 on the 
annealed survival probability, where the so-called Pascal principle will be introduced. In Section [3j 
we give a preliminary upper bound on the quenched survival probability in dimensions 1 and 2, as 
well as an upper bound for a semi-annealed system. Lastly, in Section [H we prove the existence of the 
quenched Lyapunov exponent in Theorems 11.21 and 11.31 via a shape theorem, and we show that the 
quenched Lyapunov exponent is always positive. 

2 Annealed survival probability 

In this section, we prove Theorem 11.11 We start with a proof in Section 12.11 of the existence of the 
annealed Lyapunov exponent \d,y,K,p,i>- Our proof follows the same argument as for the catalytic model 
with 7 < in Gartner and den Hollander [12], which is based on a special representation of E^[Z^] 
after integrating out the Poisson random field £, which then allows us to apply the subadditivity 
lemma. In Section 12.21 we prove Theorem 11.11 for the special case k = 0, i.e., X = 0, relying on 
exact calculations. Sections 12.31 and 12.41 prove respectively the lower and upper bound on E^[Z 7 ] 111 
Theorem 11.11 for d = 1,2 and general k > 0. The lower bound is obtained by creating a space-time 
box void of traps and forcing X to stay inside the box, while the upper bound is based on the so-called 
Pascal principle, first introduced in the physics literature by Moreau et al [191 [20]. I n Section |2~4"1 we 
will also prove the aforementioned Corollary 12.11 on the range of a symmetric random walk. 
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2.1 Existence of the annealed Lyapunov exponent 

In this section, we prove the existence of the annealed Lyapunov exponent 



A = A 



1 



lim -logE«[Z t 7 f 



(10) 



t-too t " L *'£ 

Remark. Clearly A > 0, and Theorem 1 1 . 1 1 will imply that A always equals in dimensions d = 1,2. 

For d > 3, the lower bound for the quenched survival probability in Theorem 11.21 will imply that 
A < + k < oo, while an exact calculation of A for the case k = in Section [2.21 and the Pascal 
principle in Section [2.41 will imply that A > for all 7, i>, p > and k > 0. 

Proof of (|10|) . The proof is similar to that for the catalytic model with 7 < in [12] . As in [12] , we 

can integrate out the Poisson system £ to obtain 



E*[Zy = E^[u(t, 0)] = E^E^ 



exp^-7 / £(t-8,X(s))d8 



E; 



x 



exp[uJ2(vx(t,y)-l)}], (11) 

yez d 



where conditional on X, 



v x {t,y)=E% 



exp j-7^ 6 (Y(s)-X(t-s))ds 



(12) 



with E^[-] denoting expectation with respect to a simple symmetric random walk Y with jump rate 
p and Y(0) = y. By the Feynman-Kac formula, (v x (t,y)) t>QyeZ d solves the equation 

d 

■7^v x (t,y) = pAv x {t,y) -"/5 x{t) (y)v x (t,y), y e z d t> 
«x(0,-) = 1, 

which implies that T, x (t) := X^j/ez d ^) ~ 1) ^ s * ne solution of the equation 

d 



(13) 



df 



X x (t) = - 7 v x (t,X(t)), 
Sx(0) = 0. 



Hence, T, x {t) = —7 f v x (s, X(s)) ds, and the representation (fTT]) becomes 



E^] 

We now observe that for t\,t2 > 0, 



E, 



X 



exp 



^7 / -ux(s,X(s))ds 



(14) 



(15) 



E, 



> E; 







exp< 












exp < 













J> exp < 












J> exp < 









v x (s,X(s))ds 



(s,(e tl X)(s))ds 



(16) 



= E^JE^J, 

where tl X := ((6 tl X)(s)) s >o = (X(t\ + s) — X(ti)) s >o, we used the independence of (X(s))o<. s <t 1 
and ((0t 1 X)(s))o< s <t 2 , and the fact that for s > ti, 



v x (s,X(s)) 



E 



X(s) 







exp < 








exp < 









S — tl 



5o(y(r)-A"(a-r))dr 



vo tl x(s - *i, (6» tl X)(s - h)). 



From (|16|) . we deduce that — logE^[Z^] is subadditive in t, and hence the limit in (|1U|) exists and 



^d,i,K,P,v = - sup - \ogE^[Z]A 
t>0 t 



(17) 



1. 



2.2 Special case k = 

In this section, we prove Theorem 11.11 for the case k = 0, which will be useful for lower bounding 
E^[Z^] for general n > 0, as well as for providing an upper bound on E^[Z^] by the Pascal principle. 



Proof of Theorem 11.11 for k = 0. We first treat the case 7 € (0, 00). When k = 0, (fT5|) becomes 



E € [Zjy =exp|-i/ 7 ^ u (s,0)ds 



(18) 



where Ho is the solution of (fT3j) with X = 0. It then suffices to analyze the asymptotics of uo(t, 0) as 
t — > 00. Note that the representation (fT2|) for i>o(i, 0) becomes 



(t,0) =E^[e- 7 i'o' 5o(y(s))ds ] 



(19) 



which is the Laplace transform of the local time of Y at the origin. For d = 1,2, vo(t, 0) 1 as t f 00 
by the recurrence of simple random walks, while for d > 3, t>o(t, 0) 4- Cd for some Cd > by transience. 

By Duhamel's principle (see e.g. [9l pp. 49] for a continuous-space version), we have the following 
integral representation for the solution vx of ([13p . 



(20) 



where p s (0 is the transition probability kernel of a rate 1 simple symmetric random walk on Z, d . When 
X = 0, we obtain 

vo(t, 0) = 1 - 7 / Pp.(0)«d(t - a, 0) da. (21) 
J 



y) = 1 - 7 / P P s (y - X(t - s)) v x (t - s, X(t - s)) ds, 



Denote the Laplace transforms (in t) of i>o(t, 0) and pt(0) by 



(A) = / e- M v {t,0)dt, p{\) 
Jo 



/•oo 

/ e- xt Pt (0)dt. 
Jo 



Taking Laplace transform in (|2ip and solving for vq(X) then gives 



*>o(A) 



1 



A p + jp(X/p) 



(22) 



(23) 



We can apply the local central limit theorem for continuous time simple random walks in d = 1 and 
2 (i.e., pt(0) = (^i)^ 2 (1 + o(l)) as t — > 00) to obtain the following asymptotics for p{\) as A I 0, 



p(A) 



d = 1, 



'2\ 

^1(1 + »(D), 

7T 

[ G d (0)(l + o(l)), 

with Gd{0) = fo° Pt(ty dt, which translates into the following asymptotics for £>o(A) as A 1 0: 

^27 1 



d = 2, 
d > 3, 



(24) 



fio(A) 



7 



(l + o(l)), 



d = 1, 



— •tt^tt( 1 + ^ 1 ))' d = 2 ' 

7 Aln(i) 



(25) 



1 



[p + 7G d (0) A 



(l + o(l)), 



d > 3. 
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Since vo(t, 0) is monotonically decreasing in t by (|19p . by Karamata's Tauberian theorem (see e.g. |10l 
Chap. XIII. 5, Thm. 4]), we have the following asymptotics for vo(t, 0) as t — > oo, 



'1 /2p 1 
7 V ?r 



vo(t,0) 



d = 1, 
d = 2, 
d > 3, 



(26) 



I p + 7<^(0) 

which by (|18p implies Theorem II .11 for k = and 7 € (0, 00). 
When k = and 7 = 00, we have 

E^Z^] = P(f(s,0) = Vs e [0,t]) = exp 

where ^(t, y) = P^(3 s G [0, t] : Y(a) = 0) for a jump rate p simple symmetric random walk Y starting 
from y. Note further that ip(t, y) solves the parabolic equation 

d 



at 



i>{t,y) = P Ai>(t,y), y^0,i>0, 



(27) 



with boundary conditions ip(-, 0) = 1 and ip(0, •) = 0. Therefore J2 y ez d ^(*> n) solves the equation 



— W> V) = -P A ^(i, 0) = p(l - V(i, ei)) = p<f>(t, ei), 



(28) 



where e\ = (1, 0, • • • , 0), <f>(t, e\) := 1 — ip(t, ei), and we have used the fact that X^ez d A V'(^ 2;) = 
and the symmetry of the simple symmetric random walk. Therefore 



exp < —up 



(s,ei)ds 



(29) 



By generating function calculations and Tauberian theorems (see e.g. [TTJ, Sec. 2.4] or [231 Sec. 32, 
P3]), it is known that <j)(t,ei), which is the probability that a rate 1 simple random walk starting 



from ei does not hit before time pt, has the asymptotics (j)(t,ei) = y^(l + °(1)) f° r d = 1, 
4>(t, ei ) = i^(l + o(l)) for d = 2, and 0(*,ei) = G d (0) -1 (l + o(l)) for d > 3. Therefore as t -> 00, 



logE«[Z, 7 



8pt 

Ti 

P t 



(1 + 0(1)), d=l, 



■Wn(l + o(l)), d = 2, 
lnt 



(30) 



pt 



G d (0) 



(l + o(l)), d>3, 



which proves Theorem II. II for k = and 7 = 00. 



Remark. When k = so that X = 0, the representation (|18|) allows us to easily compute the Laplace 
transform of D t := j J* * £(s, 0) ds, since E^[Z^] = E^[exp{— ftD^}]. By replacing jt with a suitable 
scale Xt/dt, where A G R, at = \ft for d = 1, = logt for d = 2, and at = 1 for d > 3, we can identify 



using the asymptotics in (|26p . As shown in Cox and Griffeath [6], applying the Gartner-Ellis theorem 
then leads to a large deviation principle for Dt with scale t/at, except that in [6], the derivation of 
Vl/(— A) was by Taylor expansion in A, which can be greatly simplified if we use the representation from 
(fT8l) instead. 



Xt 
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2.3 Lower bound on the annealed survival probability 

In this section, we prove the lower bound on E^[Z^] in Theorem 1 1.1 1 for dimensions d = 1 and 2, i.e., 
Lemma 2.1 For all 7 G (0, 00], k > 0, p > and v > 0, u>e Ziaue 



liminf ^plogE^Z, 7 ,] > -i/W— , d=l, 
liminf — logE^Z, 7 ,] > -i/tto, d = 2. 



(31) 



Proof. The basic strategy is the same as for the case of immobile traps, namely, we force the 
environment £ to create a ball B Rt of radius Rt around the origin, which remains void of traps up to 
time t, and we force the random walk X to stay inside B Rt up to time t. This leads to a lower bound 
on the survival probability that is independent of 7 G (0, 00] and k > 0. Surprisingly, in dimensions 
d = 1 and 2, this lower bound turns out to be sharp, which can be attributed to the larger fluctuation 
of the random field £ in d = 1 and 2, which makes it easier to create space-time regions void of traps. 
Note that it is clearly more costly to maintain the same space-time region void of traps than in the 
case when the traps are immobile. 

Recall that £ is the counting field of a family of independent random walks {Yj J } yeZ di < j <Ny , where 
{N y } are i.i.d. Poisson random variables with mean v. Let B r denote the ball of radius r, i.e., 
B r = {x G Z d : ||x||oo < r }- Fo r a scale function 1 « R t « \[t to be chosen later, let E t denote 
the event that N y = for all y G B Rt . Let F t denote the event that Yj'(s) ^ Bn t for all y ^ B Rt , 
1 < j < -^2/5 and s G [0, £]; furthermore, let Gf denote the event that X with X(0) = does not leave 
B Rt before time t. Then by ([3|), 

E«[Z 7 € ] > P(E t HF t n G t ) = P(E t )P(F t )F(G t ). (32) 

Note that P{E t ) = e -f(2iJt+i) . Xo estimate P(G*), note that by Donsker's invariance principle, if 
1 « Rt « y/i as t — > 00, then there exists a > such that for all f sufficiently large, 

inf P (X(s) G £^ V s G [0, t] , X(t) G B^ /2 X(0) = x) > a. (33) 
By partitioning [0, t] into intervals of length R\ and applying the Markov property, we obtain 



P(G t ) > P(X(s) G B Rt V s G [(i - l)Rf,iRt], and X(ii^) G B Rt/2 , i = 1,2,- • • , |t/fl?" 



> = e *W*?. (34) 

To estimate P(F f ), let F f denote the event that Yf(s) / for all y G Z d , 1 < j < N y , and s G [0, i]. 
Note that P(-Ft) is precisely the annealed survival probability E^[Z^] when k = and 7 = 00, which 
satisfies the asymptotics in Theorem 11.11 by our calculations in Section [2.21 We next compare P{Ft) 
with P{F t ). 

For a jump rate p simple random walk Y starting from y G Z rf , let tb r denote the stopping time 
when Y first enters B Rt , and tq the stopping time when Y first visits 0. Then standard computations 
yield 

ln¥(F t ) = -v Yl ^y( T B Rt <t), (35) 
y& d \B Rt 

and a similar identity holds for lnP(i^) with B Rt replaced by -Bo- Note that 

£ P U T B Rt <t)> £ P^(ro <*)=£ P, y (r < t) - £ F v^ < *)• 

Kt y& d \B Rt yeZ d v^B Rt 
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Hence 

lnP(F t ) < lnP(F t ) + v ^ ( r o < < InP(^) + u{2R t + l) d . (36) 

On the other hand, for e > 0, we have 

^2^(r <t + et) > F y <t,r <t + et)> Jg F y (t < et) £ P^ (r Bflt < i), 



where we used the strong Markov property. Therefore 

^ ^ ^B Rt <t)< y ~ . 

y& d \B Rt ^ou Rt z \ u _ ) 



and hence by ([35 



lnP(Ft) > (37) 



We now choose Rt for d = 1 and 2. For <i = 1, let Rt = y/t/hit, which is by no means the unique 
scale appropriate. Clearly inf z& qb F Y (tq < et) — > 1 as t — > oo. By (|3"o]) - (|57|) . the fact that P(F t ) 

yj 1 1 In £ 

satisfies the asymptotics in Theorem II .11 for k = and 7 = 00, and that e > can be made arbitrarily 
small, we obtain 

lnP(F t ) = -vJ^O- + o(l)) = lnP(F t ). 



Furthermore, for Rt = -\Jt/\n.t we have 



lnP(£ t ) = -v(2y/t/hii + 1) and lnP(G<) > In a In t, 

whence substituting these asymptotics into ()32p gives ()3ip for d = 1. 

For d = 2, let i?t = Int. Then we have inf zg aB lnt P]f (to < et) — > 1 as t — > 00, which is an easy 
consequence of |17^ Exercise 1.6.8]. By the same argument as for d = 1, we have 

ln¥(F t ) = -virp^il + o(l)) = lnP(F f ). 

Together with the asymptotics 

lnP(£ 4 ) = -i/(2 In t+ l) 2 and lnP(G t ) > 

In t 

we deduce from (|3"2"j) the desired bound in (|3"Tj) for <i = 2. I 



2.4 Upper bound on the annealed surivival probability: the Pascal principle 

In this section, we present an upper bound on the annealed survival probability, called the Pascal 
principle. 

Proposition 2.1 [Pascal principle] Let £ be the random field generated by a collection of irreducible 
symmetric random walks {Yj J }y^z d ,i<j<N y on 

Z d with jump rate p > 0. Then for all piecewise constant 
X : [0, t] — > Z, d with a finite number of discontinuities, we have 



cxp 



"7 



as,x( s ))ds 



exp 



7 / £M)ds 



(38) 
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In words, conditional on the random walk X, the annealed survival probability is maximized when 
X = 0. The discrete time version of this result was first proved by Moreau et al in [19\ I20| . where 
they named it the Pascal principle, because Pascal once asserted that all misfortune of men comes 
from the fact that he does not stay peacefully in his room. The Pascal principle together with the 
proof of Theorem 11.11 for n = in Section 12,21 imply the desired upper bound on the annealed survival 
probability in Theorem 11.11 for dimensions d = 1,2, and it also shows that for d > 3, the annealed 
Lyapunov exponent Xd. ■y,K,p,u is always bounded from below by \d,f,o,p,v — 

We present below the proof of the discrete time version of the Pascal principle from [20], which 
being written as a physics paper, can be hard for the reader to separate the rigorous arguments from 
the non-rigorous ones. We then deduce the continuous time version, Proposition 12. 1[ by discrete 
approximation. As a byproduct, we will show in Corollary 12.11 that the expected cardinality of the 
range of a continuous time symmetric random walk increases under perturbation by a deterministic 
path. 

Moreau et al considered in [20] a discrete time random walk among a Poisson field of moving traps, 
defined as follows. Let X be a discrete time mean zero random walk on Z rf with Xq = 0. Let {N y } y& %d 
be i.i.d. Poisson random variables with mean u, and let {Yj} y ^ d ,i<j<N y °e a family of independent 
symmetric random walks on % d where Yj denotes the j-th random walk starting from y at time 0. 
Let 

£(n,x):= Yl $x{Yf{n)). (39) 
y&Z d ,l<j<N y 

Fix < q < 1, which will be the trapping probability. The dynamics of X is such that X moves 
independently of the traps {X^}j,ez d ) i<j<iV w ) an d at each time n > 0, X is killed with probability 
1 — (1 — q^( n ' X ( n )\ Namely, each trap at the time-space lattice site (n,X(n)) tries independently to 
capture X with probability q. Given a realization of X, let a (n) denote the probability that X has 
survived till time n. Then analogous to (jlip . we have 

a x {n) = E« [(1 - g)E2=ofo*(0)l = eX p { - v w q ' x {n,y)\, (40) 



where if we let Y denote a random walk with the same jump kernel as Yj , then 



tt) 9,A (n, y) := 1 — E y (1 - 9 )£*=o 1 {*-(0=*(0} . (41) 

The main result we need from Moreau et al [20] is the following discrete time Pascal principle. 
Lemma 2.2 [Pascal principle in discrete time |20j] 

Let Y be an irreducible symmetric random walk on Z rf with Pq (Y(l) = 0) > 1/2. Then for all 
q £ [0, 1], n e N and X : N ->■ Z d , we have 

^2w^(n : y)> ^^' (n,y), (42) 
and hence a x (n) < a°(n), where w q '° and a denote w q ' X and a x with X = 0. 

Proof. The argument we present here is extracted from [2D]. First note that the assumption Y is 
symmetric implies that the Fourier transform f(k) := E^fe 1 ^'-^^ 1 ^] is real for all k E [— n, 7r] d . The 
assumption Pq (Y(l) = 0) > 1/2 guarantees that f(k) € [0,1]. If we let p%(y) denote the n-step 
transition probability kernel of Y, then by Fourier inversion, we have 

pl(0) > Pn(y), , .. . _ 

for all n > 0, y G Z a . (43) 

pZ(0) > p Y n+ M 

If we now regard X as a trap, then w q,x (n, y) can be interpreted as the probability that a random 
walk Y starting from y gets trapped by X by time n, where each time Y and X coincide, Y is trapped 
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by X with probability g. More precisely, let Zi, i G No, be i.i.d. Bernoulli random variables with mean 
g, where Z{ = 1 means that the trap at (i,X(i)) is open. Then X is killed at the stopping time 

T X (Y) := min{z > : Y(i) = X(i), Z, = 1}, (44) 

and w q ' X (n, y) = P^(t^ < n). 

We examine the following auxiliary quantity, where by decomposition with respect to r^-, we have 

n-l 

1 = E P ? (^ n ) = Z » = 1 )=EE P ?fo = k )ptk(X(n) - X(k)) q + P?(rjf = n) 

yez d fc=o j,ez d yez d 

71-1 

<?EE "f (rx = *)*£-*(<>) + E *f = 

k=0y£Z d y£Z d 

(45) 

where in the inequality we used (|4"3"j). Similarly, when X is replaced by X = 0, we have 

^ = E p ? = o,z n = i)= q ^2Yl *v<n = k )pt k (o) + E p ?fo = *)■ ( 46 ) 

Denote 

:= £iB**(»,y) = £ Pjfa* < n), 

y 6 Z d 2,eZ d 

j/ez d 



Note that S * = 5 ° = g, and Pf (r* = *0 = S* - S*^, E, e z d ^ffo = k) = S° k - S° k _ v where 

we set S* 1 = S?_ 1 = 0. Together with (05} and (05}, this gives 



n— 1 n— 1 



g EPn-fc(0)(5g - Sg_i) + 5° - < g E^- fc (0)(^f - + S x - S x _ v 



k=0 k=0 

Rearranging terms, we obtain 

n-2 



S x -S° n >{l- qpf (0))(S* ! - SS_i) +?E (fn-fc-i(O) -pL fc (0))(5f - Sg). (48) 

fc=0 



This sets up an induction bound for S x — S®. Since — Sq = 0, 1 — gp^(0) > 0, and Pj^(0) is 
decreasing in by (|4"3"]> . it follows that S 1 ^ > 5° for all n G No, which is precisely ([12}. ■ 

Proof of Proposition 12.11 Integrating out £ on both sides of (|38|1 as in (jlip shows that (|38p is 
equivalent to 

£>^(t,y)> ^^°(t,y), (49) 



where 



w^ x (t,y) :=1-E£ 



exp{- 7 ^\(Y( S )-X( S ))d S } 



(50) 



For n G N, let y( n )(£;) = F(f ) and X( n \k) = X(±) for it G N . Clearly yM is symmetric, and for 
n sufficiently large, Pjf™' (Y^ n \l) = 0) > 1/2. Therefore we can apply Lemma O with F = yM, 
X = X^ and g = gW = jt/n to obtain 

w^ x(n \n,y) > J2 w lt/nfi {n,y). (51) 
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By gU and the definition of and X^ n \ we have 

^/"^'"'(n y) = l-E yW ^_ 7t^£=ol {r (n) (fc)=x (n) w} = 1 _ E Y J 1 _llj^ 1 W = " u /«)} 

By the assumption that X is a random walk path which is necessarily piecewise constant with a finite 
number of discontinuities, for a.s. all realization of Y, we have 

ton fl-^ E - ol ^-^> =exp (- 7 f t S Q (Y(s)-X(s))ds\. 

n->oo V n J L Jo J 

Therefore by the bounded convergence theorem, lrm n _ > . 00 w'Y t / n > x n (n, y) = w^' X (t,y). By the same 
argument, lim n _ 5 . 00 tt) 7 '/ n,0 (n, y) = w 7 '°(t,y). Next we note that w 7 '/ n,x( ™ ) (n, y) is the probability 
that y( n ) is trapped by X^> before time n. Since and are embedded in 1" and X, we 

have w 7< / n, ^ (n) (n, y) < P y (rx < i) uniformly in n, where Tx = inf{s > : Y(s) = X(s)}. Clearly 
J2 y eZ d F y( T x < t) < oo. Similarly w^ n >°(n, y) < P^(r < t) uniformly in n and J2 ye z d F y( T < *) < 
oo. Therefore we can send n — >• oo and apply the dominated convergence theorem in (|5ip . from which 
E2) then follows. I 



The Pascal principle in Lemma 12.21 and Proposition 12.11 have the following interesting conse- 
quence for the range of a symmetric random walk, which we denote by Rt(X) = {y G Z d : X(s) = 
y for some < s < t}. 

Corollary 2.1 [Increase of expected cardinality of range under perturbation] 

Let Y and X be discrete time random walks as in Lemma \2.Sk Let Y be a continuous time irreducible 
symmetric random walk on 7L d with jump rate p > 0, and let X : [0, t] —> 7L d be piecewise constant with 
a finite number of discontinuities. Then for all n £ No, respectively t >0, we have 

E*[\R n (Y-X)\] > El[\R n (Y)\], 

E%[\R t (Y-X)\] > E%[\R t (Y)\], 

where \ ■ \ denotes the cardinality of the set. 

Proof. The first inequality in (|52p for discrete time random walks follows from the observation that 

Y, F U T x< n ) = Yl F t{Y(i)-X(i)=y for some <i < n) = E% [\R n (Y - X)\], 
y& d y& d ^ 

Y P JVo < n) = ^ P?"(y(i) = y for some < i < n) = [\Rn(Y)\] , 

y&L d y£Z d 

where t x = min{i > : Y{ = X{\ and To = min{i > : Yi = 0}, which combined with Lemma 12.21 for 
q = 1 gives precisely 

Ki^x <n)>Y < (54) 

y£Z d y&L d 

The continuous time case follows by similar considerations, where we apply Proposition 12.11 with 
7 = oo, or rather 7 > with 7 f 00. I 



3 Quenched and semi- annealed upper bounds 

In this section, we prove sub-exponential upper bounds on the quenched survival probability in di- 
mensions 1 and 2 (the exponential upper bound in dimensions 3 and higher follows trivially from the 
annealed upper bound by Jensen's inequality and Borel-Cantelli). Although they will be superseded 
later by a proof of exponential decay using sophisticated results of Kesten and Sidoravicius [16] , the 
proof we present here is relatively simple and self-contained. Along the way, we will also prove an 
upper bound (Proposition 13. 2p on the annealed survival probability of a random walk in a random 
field of traps £ with deterministic initial condition, which we call a semi-annealed bound. 
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Proposition 3.1 [Sub-exponential upper bound on Z^] There exist constants C\,Ci > de- 
pending on 7, k, p,u > such that a.s. with respect to £, we have 



logt 

i— >oo * 



lim sup ° log < — Ci , 



,. log logt 7 
lim sup log Z/ c < — O2 



d = 1, 
d = 2. 



(55) 



77ie same bounds hold if we replace by u(t,0) as in Theorem \1.3l 



Proof. The proof is based on coarse graining combined with the annealed bound in Theorem ll.il Let 
us focus on dimension d = 1 first. Let X be a random walk as in ([2]), and let M{t) := sup 0<s<t |X(s)|oo. 
The first step is to note that by basic large deviation estimates for X, 



E; 



x 



exp{ -j J Z(s,X(s))ds}l {Mt > t} \ < F*(M t > t) < e 



-ct 



for some C > depending only on k. Therefore to show ([55]) . it suffices to prove that 



E 



-V 



exp{- 7 ^ £(s,X( S ))ds}l 



{M t <t} 



Ct 

< e lo s* 



(56) 



for some C > for all t sufficiently large. Since the integrand in the definition of is monotone in 
t, we may even restrict our attention to t E N. 

The second step is to introduce a coarse graining scale Lt := A logt for some A > 0, and partition 
the space-time region [— 2t, 2t] x [0, t] into blocks of the form Aj & := [(i — l)L t , iL t ) x \{k — l)Lf, kh\) 
for i,k eZ with -j^ + l<i<j^ and 1 < k < j?. We say a block 



K i>k is good if Z((k - l)L 2 t ,x) > 

(i—l)Lt<x<iLt 



Since for each s > 0, (£(s,x)) x& z are i.i.d. Poisson distributed with mean v, by basic large deviation 
estimates for Poisson random variables, there exists C > such that for all t > 1, 

P(A i>k is bad) < e~ CuLt . 

Let Gt(£) be the event that all the blocks Aj & in [— 2t, 2t] x [0,t] are good. Then 

which is summable in t > 2, t E N, if j4 is chosen sufficiently large. Therefore by Borel-Cantelli, a.s. 
with respect to £, for all t E N sufficiently large, the event Gt(£) occurs. To prove ()55|) . it then suffices 
to prove 



l Gt (^o [exp { - 7 jf £(s, ds} \ {Mt 



<*} 



< e lo s* 



(57) 



almost surely for all i E N sufficiently large. 

The third step is applying an annealing bound. More precisely, to show (|57p . it suffices to average 
over £ and show that 



E^Ejf 



expj-7^ £(s,X(s))ds}l {Mt<t} l Gt 



(0 



2Ct 

< e lo e* 



(58) 



for some C > for all t E N sufficiently large. Indeed, (|58|) implies that 



Pf (l Gt ( E? [exp { - 7 y* £( S , *(«)) ds} l {Mt 



<*} 



Ct 



Ct 



> e lo s * < e lo s * . 
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from which (|57p then follows by Borel-Cantelli. 

To prove (|58|) . let us denote := exp { — 7 f^ k J^ L 2 C( s i A(s)) ds}, and let be the cr-field 
generated by (X s ,£(s, -))o< s <kL 2 - Replacing by t/[t/L 2 \ if necessary, we may assume without loss 
of generality that t/L 2 t = t/{A\ogt) 2 G N. Then 



E^Eff 



E^Eff 



exp I -1 J C(s,X(s))ds^l{M t <t}^G t 



(0 



E S E, 



A 



t/L? 
fc=l 



i{M t <t}iGt(o n 



fe=i 



(59) 



By Proposition 13.21 below, on the event \X{{k — )|oo < t and A^fc is good for all — jj- + 1 < i < 
which is an event in J- k -i, we have 



E*E;*[Z fc |.F fe _ 1 ] = E«E* 



exp 



f e( a ,x( a ))d S } 

(fe-l)Lf J 



< e 



-CL t 



(60) 



for some C > depending on 7, k, p, za Substituting this bound into (j59|) for 1 < < i/L? and using 



the fact that Ilffi 



is a martingale then gives the desired bound e Ct / Lt = e 



•Ct/(A logt) 



E€E*[Z fe |.F fe _i] 
for (p|). 

For dimension d = 2, the proof is similar. We choose L( = ^4 logt with A sufficiently large. We 
partition the space-time region [— 2t, 2t] 2 x [0, t] into blocks of the form Aj^^ := [(i — l)L t ,iL t ) x [(j — 
l)L t ,jL t ) x [(fc— fcij), and we define good blocks and bad blocks as before. Applying Proposition 
32] below then gives an upper bound of exp { — Cj? ,^ } = exp{ — log A ^ og i ogt }, analogous to ([58]). 

Lastly we note that the arguments also apply to the solution of the parabolic Anderson model 



(t, 0) = E* [ exp { - 7 J £(t - s, X(s)) ds} 



The only difference lies in passing the result (|55p from t G N to t G K, due to the lack of monotonicity 
of u(t, 0) in t. This can be easily overcome by the observation that for n — 1 < t < n with n G N, 

and the fact that almost surely £(r, 0) dr < \fi for all i large by Borel-Cantelli because £(r, 0) dr 
has finite exponential moments. I 



The following is a partial analogue of Theorem ll.il for £ with deterministic initial conditions. 

Proposition 3.2 [Semi-annealed upper bound] Let £ 6e defined as in (CJ) deterministic 
initial condition (£(0,x)) xe %d. For L > and i = (ii,--- , i^) G Z d , let B L ^ := [(i\ — l)L,i\L) x 
• • • x \{i-d ~ 1-)L,idLi). Assume that there exist a > 2 and v > suc/i t/iat /or all i G [— 3L a , 3L a ] d , 
Sx-pr -£(0> x ) — Then there exist constants Cd > 0, d > 1, such that for all L sufficiently large 

and for all x G Z d with \x\oo < L, we have 



E*Ef [exp {-7^ £( S ,A( S ))ds} 



< ^ 



-CiL 



e C2 i««i. 



d = 1, 

d = 2, 
d > 3. 



(61) 



The same is true if we replace Jq £(s,X(s))ds by ^(s,X(L 2 — s))ds. 
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Proof. The basic strategy is to dominate (£(L 2 /2, 2;))| a; | 00 <2L a from below by i.i.d. Poisson random 
variables, which then allows us to apply Theorem 1 1.11 We proceed as follows. 

Let £ be generated by independent random walks {Yj) y ^z d ,x<j<(,{0,y) as m ©> an d let £ be gen- 
erated by a separate system of independent random walks (Yj) yeZ d i<j<£(o,y)i where (£(0, y)) y£ z d are 
i.i.d. Poisson distributed with mean v. Choose any v G (0, v). Then by large deviation estimates for 
Poisson random variables, 



^{G\ 



pi 



< 



E e"(o, : 
pi 



E 

i£{-3L a ,3L a ] d 



> for some ^ [-3L a ,3L a ] 

£(0,x) > vlf 



(62) 



E 

x6-B r 



l<i<?(0,J/) anC ^ (yf)y£l d ,l<j<£(0,y) aS 

follows. For each walk E/ with 1 < j < £(0, y) and y G i? L -? for some z G [— 3L a , 3L a ] d , we can match 

: k iui Duinc & c ^ L j- and 1 < k < £(0, z), which is possible on the event Gl- 



On the event Gl, we will construct a coupling between (E 7 - ) 

E^ with a distinct walk Ejf for some z £ B L 

Independently for each pair of walks (E^,E fc 2 ), we will couple their coordinates as follows: For 
1 < i < d, the i-th coordinates of the two walks evolve independently until the first time that their 
difference is of even parity. Note that this is the case either at time already or at the first time when 
one of the coordinates changes. From then on the i-th coordinates are coupled in such a way that 
they always jump at the same time and their jumps are always opposites until the first time when 
the two coordinates coincide. From that time onward the two coordinates always perform the same 
jumps at the same time. For walks in the £ and £ system which have not been paired up, we let them 
evolve independently. Note that such a coupling preserves the law of £ (resp. £), and each coupled 
pair (Yj, E fe 2 ) is successfully coupled in the sense that Yj / (L 2 /2) = Y£ {L 2 /2) if the trajectory of Yj is 
in the event 



{ sup (E/(i)-Ef(0))i G [L/2,L] and inf (Y*(t)-Y?(0)) i G [-L, -L/2] V 1 < i < d] 

U^f2/o J J 0<t<L 2 /2 J J i 



0<t<L 2 /2 



because \y — < L by our choice of pairing of Yj and E fc 2 . Then by our coupling of £ and £, on the 
event Gl, we have 



£(L 2 /2,x) > Q(x) := J2 



{yy(Ly2)=x} 



for all Ixloo < 2L a . 



(63) 



i<i<£(o,») 



Now observe that, because (£(0, x)) xeZ d are i.i.d. Poisson with mean u, and (Xf) y eZ d ,l<j<^(0,y) are 
independent, (C( x )) X £Z d are a ^ so i-i-d- Poisson distributed with mean a := v Ft(Ef) = i/P^-E?), which 
is bounded away from uniformly in L by the properties of simple symmetric random walks. This 
achieves the desired stochastic domination of £ at time L 2 /2. Let Cl(*)") denote the counting field 
of independent random walks as in ([1]) with initial condition Cl(0, y) = Qiy)^{\ y \ ao <2L a }- Then using 
uniformly in x G Z rf with |x|oo < L, we have 



E € E 



A 



exp{-7jf £ e(s,X( S ))d s } 



E 5 '% 



x 



exp 



< pZ(G c L ) + FX(\X(L 2 /2 



>L 2 



+ sup E Ci Ef 

Moo<L 2 



< 



6 d L a de -C^L d +e -CL* + 



sup E^E* 

\x\ x <L 2 



Hit 

exp | — 7 

■L 2 /2 



t(s,X(s))ds 

L 2 /2 



( L (s,X(s))ds} 



exp 



{"7^ C L (8,X( 8 ))d8}] 



By the same argument as for (jlip . we have 

r-L 2 /2 



exp{- 7 / CL(s,X(s))ds}] =Ef [exp{-a ^ (1 - t;x(£ 2 /2, y))} 

•^0 I.J xnrn 



(64) 



(65) 



Moo<2L° 
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where 



v x (L 2 /2,y) =E 



Y 



exp 



7 



L 2 /2 



,5o(y( s )-x( s ))d S } 



To bound ()65|) . note that by a union bound in combination with Azuma's inequality we obtain, 



sup P*( sup |X(s)|oo > 2L 2 ) < e" 

Moo<L 2 0<s<L 2 /2 



OL 2 



(66) 



On the complementary event {sup < s <£2/ 2 |^(s)|oo < 2L 2 }, we have 

1 - v x (L 2 /2,y) < P^(r 2L2 < L 2 /2) < P(P L 2 /2 > |y| 



2L^ 



where r 2 ^2 := inf{s > : |Y(s)|oo < 2L 2 }, and Vip.^ is a Poisson random variable with mean pL 2 /2, 
which counts the number of jumps of Y before time L 2 /2. Therefore for L sufficiently large, 

(l-v x (L 2 /2,y))< ]T ¥(V L 2 /2 >\y\ O0 -2L 2 ) 

\y\ aa >2L" \y\ oc >2L a 



< C nVL V 2>rl2)r d - 1 < CM[P k L y 2 ] ^ r^ 1 

r=2L a 

< C(pL 2 /2) k (2L a 



r=2L a r=2L a 
- a\d—k ^ £j j^—(a—2)k+ad ^ ^ 



(67) 



where we have changed the values of the constant C (independent of L) from line to line, and the last 
inequality holds for all L large if we choose k large enough. Substituting the bounds (|66|) - (|67|) into 
(|65|) then gives the following bound uniformly for x G Z rf with |x|oo < L 2 : 



exp 



L 2 /2 







( L (s,X( S ))ds} 



< 



sup \X(s)\ OQ >2L 2 ) +E* expl -aY (l-v x (L 2 /2,y)) + l 

0<s<L 2 /2 1 1 



< e 



exp { - a^2(l ~ vx(L 2 /2,y))}], 

y&L d 



where by the representation (|lip . the expectation is precisely the annealed survival probability of a 
random walk among a Poisson field of traps with density a, for which the bounds in apply with v 
replaced by a and t by L 2 /2. Substituting this bound back into (|64p then gives (|6ip . The same proof 
applies when we reverse the time direction of X in (|6ip . I 



4 Existence and positivity of the quenched Lyapunov exponent 

In this section, we prove Theorems 11.21 and 11.31 In Section 14.11 we state a shape theorem which 
implies the existence of the quenched Lyapunov exponent for the quenched survival probability Z^. 
In Section 14.21 we prove the stated shape theorem for bounded ergodic random fields. In Section 14.31 
we show how to deduce the existence of the quenched Lyapunov exponent for the solution of the 
parabolic Anderson model from what we already know for the quenched survival probability. Lastly 
in Section 231 we prove the positivity of the quenched Lyapunov exponent, which concludes the proof 
of Theorems 11.21 and 11.31 

4.1 Shape theorem and the quenched Lyapunov exponent 

In this section, we focus exclusively on the quenched survival probability Z^. The approach we adopt 
in proving the existence of the quenched Lyapunov exponent for Z^ uses the subadditive ergodic 
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theorem and follows ideas used by Varadhan in [25] to prove the quenched large deviation principle 
for random walks in random environments. 

For s > and x G Z d , let P^f s and E^ s denote respectively probability and expectation for a jump 
rate k simple symmetric random walk X, starting from x at time s. For each < s < t and x, y £ Z d , 
define 



e(s,t,x,y,g) := E: 



f 

x 

x.s 



exp <^ -7 / X{u)) du \ l {X (t)= y } 



(68) 



a(s, t, x, y, f) := - log e(s, t, x, y, £). 

We call a(s,t,x,y,£) the point to point passage function from x to y between times s and t. We will 
prove the following shape theorem for a(0, t, 0, y, £). 

Theorem 4.1 [Shape theorem] There exists a deterministic convex function a : M d — > M, which 
we call the shape function, such that W^-a.s., for any compact K C M. d , 

lim sup |t _1 a(0,t,0,y,0 - a(y/t)| = 0. (69) 
t_> °° ydtKr\Z d 

Furthermore, for any M > 0, we can /ind a compact K cM. d such that 

1 



lim sup - log E(f 

t— >-oo t 



exp | -7^ f(s,X(s))dsl l{x(t)^ti<r} 



< — M. (70) 



Remark. Note that ([5]) in Theorem 11.21 follows easily from Theorem 14. 11 which we leave to the reader 
as an exercise. In particular, the quenched Lyapunov exponent satisfies 

Vy W = inf a{y) = a(0) = lim r^O, t, 0, 0, 0, (71) 

where inf^gjgd a{y) = a(0) follows from the convexity and symmetry of a, since £ is symmetric. 

The unboundedness of the random field £ creates complications for the proof of Theorem 14.11 
Therefore we first replace £ by £at := (max{£(s, x), N}) s>0 ^ x£Z d for some large N > and prove a cor- 
responding shape theorem, then use almost sure properties of £ established by Kesten and Sidoravicius 
in |15j to control the error caused by the truncation. 



Theorem 4.2 [Shape theorem for bounded ergodic potentials] Let ( := (C( s > x )) s >o xez d be a 
real-valued random field which is ergodic with respect to the shift map 9 r ,zC := (C( s + r > x + z )) s >o zez d > 
for all r > and z 6 7L d . Assume further that |C(0, 0)| < A a.s. for some A > 0. T/ien i/ie conclusions 
of Theorem\4-l\ hold with £ replaced by £. 



Remark. Note that Theorem 14 . 2 1 can be applied to the occupation field of the exclusion process or the 
voter model in an ergodic equilibrium, which in particular implies the existence of the corresponding 
quenched Lyapunov exponents. 

Before we prove Theorem 14.21 in the next section, let us first show how to deduce Theorem 14.11 
from Theorem 14. 2\ using almost sure bounds on £ from [T5] . 

Proof of Theorem 14.11 Note that, since £ is non- negative, ()70|) follows from elementary large 
deviation estimates for the random walk X, if we take K to be a large enough closed ball centered at 
the origin, which we fix for the rest of the proof. 

By applying Theorem 14.21 to the truncated random field we have that for each TV > 0, there 
exists a convex shape function : M d — )• M such that ([69]) holds with £ replaced by £tv and a replaced 
by «tv. Note that cxn is monotonically increasing in N, and its limit a is necessarily convex. To prove 
(|69p . it then suffices to show that, for any e > 0, we can choose N sufficiently large such that P^-a.s., 

- sup |a(0, t, 0, y, £) — a(0, t, 0, y, £n)\ < e for alH sufficiently large. (72) 

* y£tKnZ d 
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To prove (|72[) . we will need Lemma 15 from |15] . which by Borel-Cantelli implies that there exist 
positive constants Co, Ci, C2, C3, C4 with Co > 1, such that if Ej denotes the space of all possible 
random walk trajectories 7r : [0, t] — > Z d , which contain exactly I jumps and are contained in the 
rectangle [— C\tlogt, C\tlogi\ d , then P^-a.s., for all t G N sufficiently large, we have 



sup / £(s,ir(s))l 
vr£S ; Jo 



vC d m} ds < (t + jr C 3 Co (d+6)+d e- C4C o /4 Vm£N,l>0, (73) 



where y4 m := C 3 C ,(d+6)+ V^ /4 as m -> 00. 

One important consequence of ([75]) is that 

< sup a(y) < 00. 

Indeed, if h(X) denotes the number of jumps of X on the time interval [Q,t], then 



(74) 



sup a(y) < lim —t 1 log inf E 



x 



y£K 



t— >co 



y&Knl 



X 



< lim —t 1 log inf E 

t->oo y&tKr\Z d 







exp < 








exp < 









7 / f («, *(<0) ds f l{X(t)=y,« t (X)<2D(A-)t 



where D(K) := sup^g^- |y|i. We can then apply ([73]) and large deviation estimates for random walks 
to the above bound to deduce sup y£K a(y) < 00. The fact that sup y£K a(y) > for a large ball K 
again follows from basic large deviation estimates. 

By large deviation estimates, we can find B large enough such that 

^o(lt{X) > Bt or X £ 3 k ( X )) < e - 2sup «£* Q ^* for all t sufficiently large. 

Let N = C 2 vC^ m . Then by {73]), P^-a.s., uniformly in y G Z d and for all t large, we have 



(75) 



e(0,i,0,y,£) > e -(i+^7* E £ exp |_ 7 J* ^( S) X (s)) dsX 1 
> e-( 1+B )^(e(0,t,0,j/,^)- 



{X(t)=y,lt(X)<Bt,X&S lt{x) } 



-2 sup 



a(y)t 



(76) 

where in the last inequality we applied ([75]) . Since — i -1 log e(i, 0, y, £tv) — > a^{y/t) uniformly for 
y G tKnZ d by Theorem 14. 2| and sup ygA - a^{y) < sup y£K a(y), ([76]) implies that P^-a.s., uniformly 
in y G tK n Z rf and for all £ large, we have 

rVo, t, 0, y, e) < i _1 a(0, i, 0, y, £ N ) + (1 + B)A m7 + o(l). 

Since a(0, £, 0, y, £) > a(0, t, 0, y, £jv), and A m can be made arbitrarily small by choosing m sufficiently 
large, ([72]) then follows. I 



Remark. Theorem 14.11 in fact holds for the catalytic case as well, where we take 7 < in (|68|) and 
(|70p . This implies the existence of the quenched Lyapunov exponent in Theorem 11.21 for the catalytic 
case, where we set 7 < in the definition of Zj^. Indeed, Theorem 14.21 still applies to the truncated 
field £at. To control the error caused by the truncation, the following modifications are needed in the 
proof of Theorem 14.11 To prove ([70]) , we need to apply ([73]) . More precisely, we need to first consider 
trajectories (X x )o< s <t which are not contained in [— C\t log t, C±t log t] d . The contribution from these 
trajectories can be shown to decay super-exponentially in t by large deviation estimates and a bound 
on £ given in (2.37) of Lemma 4]. For X which lies inside [—C\t\ogt, C\t log t] d , we can then use 
([75]) and large deviations to deduce ([TP]) . In contrast to ([75]) . we need to upper bound e(0, t, 0, y, £) 
in terms of e(0, t, 0, y, £jy). The proof is essentially the same, except that in place of ([75]) . we need to 
show that we can choose B large enough, such that P^-a.s., 



sup 



exp ^ |7| / £(s, X(s)) ds f l{x(t)=w} 



)>Bt or x^- 



< inf F*(X t = y). (77) 



This can be proved by appealing to ([TO]) , and applying ([75]) and large deviation estimates. 
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4.2 Proof of shape theorem for bounded ergodic potentials 

In this section, we prove Theorem 14.21 From now on, let Q+ denote the set of positive rationals, and 
let Q d denote the set of points in M. d with rational coordinates. We start with the following auxiliary 
result. 

Lemma 4.1 There exists a deterministic function a : Q d — > [—7^,00) such that for every y G Q d , 

lim t _1 a(0,t,0,ty,C) = a(y) P f - a.s. (78) 

t— >oo 
tySZ d 

Proof. Since we assume y G Q d and ty G Z d in (|78p . without loss of generality, it suffices to consider 
y G Z d and t G N. Note that by the definition of the passage function a in (|68p . P^-a.s., 

a(t 1 ,t 3 ,x 1 ,x 3 ,() < a(t 1 ,t 2 ,x 1 ,x 2 ,C) + a(t 2 , t 3 , x 2 , x 3 , Q Vti < t 2 < t 3 , Xi,x 2 ,x 3 G Z d . (79) 

Together with our assumption on the ergodicity of C, this implies that the two-parameter family 
a(s,t, sy,ty,(), < s < t with s,t 6 Z, satisfies the conditions of Kingman's subadditive ergodic 
theorem (see e.g. HE]). Therefore, there exists a deterministic constant a(y) such that (fT8"|) holds. 
The fact that a(y) > — "fA follows by bounding £ from above by the uniform bound A, and a(y) < 00 
follows from large deviation bounds for the random walk X. I 

To extend the definition of a(y) in Lemma |4. II to y ^ Q d and to prove the uniform convergence in 
we need to establish equicontinuity of t~ 1 a(0, t, 0, ty, () in y, as t — > 00. For that, we first need 
a large deviation estimate for the random walk X. 

Lemma 4.2 Let X be a jump rate k simple symmetric random walk on Z d with X(0) = 0. Then for 
every t > and x G Z rf , we have 

— J(-) t 

F*(X(t) =x)= d _* ± /4 (1 + o(l)) , (80) 



(2 7 rt ) fnt 1 (^ + S 
where 

J ( x ) '■= ^ with j(y) := ysinh" 1 ?/ - ^ y 2 + 1 + 1, 

i=i 

and the error term o(l) tends to zero as t — >• 00 uniformly in x G i-fT H Z d ; /or any compact K C M d . 

Proof. Since the coordinates of X are independent, it suffices to consider the case X is a rate re/d 
simple symmetric random walk on Z. Let a := . Let Z\, • • • , Zi^-i be i.i.d. with 

P(Z A = y)= F{X(a) = y)e^*W, y G Z, 

where 

$(A) = logE[e AX ( CT )] = ^(coshA - 1). 

a 

Note that 

E[Z?] = ^?(A) = ^ sinhA and Var(Z 1 A ) = ^|(A) = ^ cosh A. 

We shall set A = sinh _1 (ff ) so that E[Z X A ] =x/\t\. If we let S m := Yl\=i Z?, then observe that 
F$(X(t) = x) = ¥(S m = x ) e -^+m*(A) = P (5 m = x ) e -5i(^)*. 



/ 2 2 

Note that — x has mean 0, variance t\i + and characteristic function 



\t~] ($(ife+A)-$(A)) -iifcx _ ia:(sinfc-fc)-ty^ + ^(l-cosfe) 

Applying Fourier inversion then gives (|80|) . 

With the help of Lemma 14.21 we can control the modulus of continuity of t~ 1 a(0, t, 0, ty, £). 
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Lemma 4.3 Let K be any compact subset of Mr. There exists 4>k ■ (0, oo) — > (0, oo) with lim^o <^x( r ) 
0, such that for any e > 0, F^-a.s., we have 



limsup sup t 1 |a(0, t, 0, x, () — a(0, t, 0, y, ()\ < 4>k((-)- 

t— >CO x,y£tKnZ d 

\\x—y\\<£t 



(81) 



Proof. Let K C M. d be compact. It suffices to consider e G (0, 1/2), which we also fix from now on. 
First note that, by Lemma 14.21 P^-a.s., 



inf e(0, t, 0, z, () > e~ m inf P£ (X t = z) > e 



-(A+l)t-sup ueK J(u)t 



for all t sufficiently large. 

Also note that for all z E Z d and t > 0, 



e(0, t, 0, z, C) = e(0, (1 - e)t, 0, C) e((l - e)t, t, w, z, 0- 



(82) 



(83) 



By large deviation estimates, we can choose a ball Br centered at the origin with radius R large 
enough and independent of e, such that K C Br and P^-a.s., 



sup ^2 e (°' ( X - e)£,0,w,C)e((l - e)t,t,w,z,() 



< P* (X((l - e)t) £ tJ3 B ) e A< < e -(^+2)i-su PtieA . * 

for all t sufficiently large. In view of (|82p . the dominant contribution in (|83p comes from u; G ti?R nZ d . 
Therefore to prove (|8ip . it suffices to verify 



limsup sup i 1 

t— >co x, y etKnx d 



log 



T, w& B R md e (°> C 1 - e K °> w > e((l - e)t, i, to, y, C) 



E^et^n^ e (°> (! - °> w ' C) e((l - e)t, t, w, x, C) 

Note that P^-a.s., and uniformly in x, y G tiC n Z rf with ||x — < et, 

EwetBanz" e(0» C 1 ~ e K 0. w > e((l - e)t, t, w, y, Q 
l>2w£tB R nz d e (°> i 1 ~ e )*> °> w > C) e((l - e)t, t, w, x, C) 

e((l-e)t,t,w,y,Q , {X (et) = y - w) 2Aet 

w&B R nZd e((l - e)i, t, w, x, C) ™ etBi? nzd ( x ( et ) = z - w) 



< ^(e). (84) 



( (x — w 

< exp < et sup I J I 

wetB R nz d ^ ^ et 



J 



ct 



+ 3Aet 



< exp^et sup | J(u) - J(v)\ + 3Aet > 

[ U,v£B 2R / e ,\\u— v\\<l J 

for all t sufficiently large, where we applied Lemma I3~2l and B 2 r/ 6 denotes the ball of radius 2R/e, 
centered at the origin. Therefore (|84p holds with 

(j) K (e) = 3Ae + e sup | J(u) - J(v)\. 

u,v£B 2R/€ ,\\u-v\\<l 



It only remains to verify that 4>k{c) i as e \. 0, which is easy to check from the definition of J. 



Proof of Theorem 14.21 Because £ is uniformly bounded, (|7U|) follows by large deviation estimates 
for the number of jumps of X up to time t. Lemma 14.31 implies that for each compact K C M. d , the 
function a in Lemma 14. 11 satisfies 

sup \a{u) — a(v)\ < 4>k(c) for all e > 0. (85) 

u,v£Kn(} d 
\\u-v\\<e 
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This allows us to extend a to a continuous function on R . 
To prove ([69]) . it suffices to show that for each 5 > 0, 

limsup sup |i _1 a(0, t, 0, y, £) — a(y/t)\ < S. (86) 

We can choose an e such that 4>k(^) < We can then find a finite number of points x±, ■ ■ • , x m G Q d 
which form an 6-HGt in li. j and along ct subsequence of times of the form t n = na with oxi G Z d for all 
Xj, we have i~ 1 a(0, f n , 0, t n Xj) — > a(xi) a.s. The uniform control of modulus of continuity provided by 
Lemma 14.31 and (|85p then implies (|86p along i n . This can be transferred to t — > oo along R using 

e(0, t, 0, y, C) > e(0, a, 0, y, () e(s, t, y, y, Q > e(0, s, 0, y, C) e -^ A ^-^ for s < *. 

Lastly, to prove the convexity of a, let x,y G R rf and /3 G (0, 1). Then P^-a.s., we have 

a(0, t n ,0,Py n + (1 - /3)x n ,C) < a(0, /3t n ,0, /3y n ,C) + a(/3t n ,t n , f3y n , (3y n + (1 - P)x n ,(), 

where we take sequences i n , x n ,y n with t n — >■ oo, x n /t n — > x, y n /t n — > y, and fly n , (1 — /3)x n G By 
Lemma |4"TT| the first term divided by t n converges a.s. to a((3y + (1 — f3)x), the second term divided 
by (3t n converges a.s. to a(y), while the last term divided by (1 — f3)t n converges in probability to a(x) 
by translation invariance. The convexity of a then follows. I 



4.3 Existence of the quenched Lyapunov exponent for the PAM 

Proof of (|9D in Theorem 11.31 Since is equally distributed with u(t, 0) for each t > 0, 

converges in probability to the quenched Lyapunov exponent \d,i,n,p,v- It only re- 
mains to verify the almost sure convergence. We will bound the variance of logu(t,0), which is the 
same as that of logZ^, and then apply Borel-Cantelli. 

Assume that t G N. Note that we can write £ as a sum of i.i.d. random fields (£i(s,x)) s>Q xeZ d, 
each of which is defined from a Poisson system of independent random walks with density u/t, in the 
same way as £. Then we can perform a martingale decomposition and write 

t t 

log zr tA - [log zr tA \ = £ y -=£ (^[iog z^ih, •••,&]- e s [log z^ia, • • • , e, 



8=1 



i=l 



i-1 



and hence Var(logZ^) = £- =1 E«[lf]. 
For each 1 < i < t, we have 

Vt = E^'-Aflog^-E^pog^,]] 



E £i+i,-,6 E £i 



log : 



E x ^-7/0 (Ei< i < t ,^f J (s,x( S ))+ei( S ,x( s ))) d,j 



= E &+i,-,& E £ Hog E ^[ e -7 Jo* ds ] _ logE^fe-T^fK^^W)^] 

where £^ denotes an independent copy of £j, and E X '* denotes expectation with respect to the Gibbs 
transform of the random walk path measure Pq', with Gibbs weight e _7 ^° ^i<j<'j^^ s,A 'W ds . Then 
by Jensen's inequality, 

E«[tf] < E^[(logE^[e- 7 /o^(^(«))^] _ logE X,i[ e - 7 ti? i (s,X(s))ds]f 

T,x,i r P -7/oe-( s ^( s )) ds ]) 2 



< 2E^ (logE^e-^^OO)^]) 2 +2 E^ (logE^fe 



E 



A'. 



= 4E« 
< 4E C 



< 47 2 E ? E Xli 



logE^fe-TJIffi^W)^]) 2 ' 



7 / &(s,X(s))d S 



4 7 2 E X '*E^ 



t e 4 (s,x( s ))d S N ' 



o 
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where in the third line we used the exchangeability of and in the fourth line we applied Jensen's 

inequalitjQ to the non-negative convex function — logx on the interval (0, 1]. 
Note that for any realization of ((X(s))o< s <t, we have 

E 6 [(/*&(s,X(s))d S ) 2 ] =2 J! E?*[£i( U ,X{u))ti(v,X(v))]dudv 

0<u<v<t 



V 2 _ v ,V 2 Vs 



£ ¥ lu(Y(v) = X(v)) + + j)^ )>u {Y(v) = X(v))) dudv 



o< u <v<t »ez; , 

ft 



< 2v 2 + 2vl P^ (y(s) = 0)ds, 
Jo 

where P^ s denotes probability for a simple symmetric random walk on Z rf with jump rate p, starting 
from y at time s, and in the last line we used that P^ (Y(s) = y) is maximized at y = for all s > 0. 
Combined with the previous bounds, we obtain 

* /"* 
Var(logu(t,0)) = Var(logZ t 7 £ ) = S^E^V 2 } < 8 7 Vi + 8 7 2 z^ / F% (Y(s) = 0) ds < Ct* 

i=i J o 

for some C > 0, since JqFq (Y(s) = 0) ds is of order \Jt in dimension d = 1, of order logt in d = 2, 
and converges in d > 3. Therefore for any e > 0, 

C 

logu(t,0) -E s [logu(t,0)]| > et) < 



which by Borel-Cantelli implies that along the sequence t n = n 3 , n G N, we have almost sure conver- 
gence of — t~ l log u(t, 0) to the quenched Lyapunov exponent \d,y,K,p,w 

To extend the almost sure convergence to t — >■ oo along R, consider i G [t n ,t n+ i) for some n G N. 
As at the end of the proof of Proposition 13.11 we have 

u(t,0) > e- K(i -*" ) e- 7 ^« (s ' 0)ds n(t„,0), 
n(t,0) < e K ( t "+ 1 -*)e 7 ^ ,I+1 ^ s ' 0)ds n(t n+ i,0). 

Note that (t n +i — t n )/t n — > as n — > oo, and we claim that also t^ 1 J t * n+1 £(s, 0) ds — >• a.s. as n — > oo, 
which then implies the desired almost sure convergence of t _1 logu(t,0) as t — > oo along R. Indeed, 
since £(s, 0) ds has finite exponential moments, as can be seen from (|15p applied to the case 7 < 
and X = 0, we have exponential tail bounds on J £ (s, 0) ds, which by Borel-Cantelli implies that a.s. 
sup <j <m £( s ' 0) < log m for all m G N sufficiently large. The above claim then follows. I 



4.4 Positivity of the quenched Lyapunov exponent 

In this section, we conclude the proof of Theorems 1 1 . 2 1 and 1 1 . 3 1 by showing that the quenched Lyapunov 
exponent Xd,y,K,p,u is positive in all dimensions. The strategy is as follows: Employing a result of 
Kesten and Sidoravicius |16} Prop. 8], we deduce that P^-a.s. for eventually all integer time points 
t, sufficiently many X paths encounter a ^-particle close-by for of order t many integer time points. 
Using the Markov property, we then show that with positive Pq' probability, X moves to a close-by 
^-particle (which itself stays at its site for some time) within a very short time interval and collects 
some local time with this ^-particle. This then implies the desired exponential decay. 

Proof of Theorems 11.21 and 11.31 Since we have shown the quenched Lyapunov exponent \d,^,K,p,u 
in Theorems II .21 and II. 31 to be the same, it suffices to consider only Theorem 11.21 Note that the upper 

1 Note that this is where the proof fails for the 7 < case. 
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bound on Xd.-y,K,p,u m Theorem 11.21 follows trivially by requiring the walk X to stay at the origin. To 
show Xd,j,K,p,v > 0] we wm make the strategy outlined above precise. In compliance with [16] we let 
Co and r > be large integers and for i G Z d define the cubes 

d 

3=1 

In a slight abuse of common notation, let -D([0, oo), Z d ) denote the Skorohod space restricted to those 
functions that start in at time and have nearest neighbour jumps only. Then set 

Jk ■= G D([0, oo),Z d ) : $ jumps at most cICq(k V l)k times up to time k}. 

For integer times t > define 

5(t):= fl J *" 

k=[t/4\ 

Then standard large deviation bounds yield 

F$(X G ~{t) c ) < e - c ('+°W), (87) 
for some c > 0. In addition, define the cube 

Ct ■= [-dc^(KVi)t,dC^(KVi)t] d nz d , 

as well as for arbitrary t G N, k G {0, • • • , t}, $ G H(i) and e > the events 

A(t,$,fc,e) := {3?G C t : $(fc) G Q r (i*) and 3y G Q r (?) : f(s,y) > lVs G [fc.fc + e/p]} 

and 

:= n { e ^(t^M > 

*6H(t) fce{|f/4j,-,t-l} 

which both depend on £. 

For e small enough, using Borel-Cantelli, it is a consequence of \16\ Prop. 8] that P^-a.s., G(t) occurs 
for eventually all t G N. Indeed, denoting by n t / 4 i ... t \ the subset of (Z d ){L*/ 4 J' - "'*} obtained by 
restricting each element of to the domain { [_i / 4J , • • • , t}, we estimate 

P«(G(t) c )<P«( |J { ^ l Wl£ )<4) 
*€H(t) fee{Lt/4j,-,t-i} 

<^( U { E Wo)<5}) 

*es(t) fee{LV4j,-,t-i} 

+ I S (*)I{L*/4J,-,*}I x ™-^ pe ( E U(t,$,fc,o)^> E U(t,$,M^ e *) 
~ u fce{L*/4J,-,t-i} fce{L*/4j ,-,*-!} 

< P? (U{ E lA( W .,o)<2}) + l H ( t )l{LV4j,-,i}l xP (E^^ et )' ( 88 ) 
*6H(t) fee{|t/4j,-,t-i} j=l 

where in the last step we observed that, given $ G by the strong Markov property of £ applied 
successively to the stopping times n := inf{j > [t/4:\ : Xlfc=Li/4j l-A(t,$,fc,o) = *}> we can couple £ with 
a sequence of i.i.d. Bernoulli random variables (pi j£ )j e N with 

P(pi, e = 1) = P^y^s) = OVs G [0,e/p] | £(0,0) > l), 
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t/2 

such that l A (t,<s>, n ,e) > Pi,e a.s. for all i G N, and hence Efce{|t/4],-,t-i} W,*,M > Ei=iPi,e on tne 
event Efce{L*/ 4 J - t-i} ^-A(t,^,k,o) — \- Here pi je corresponds to the event that given A(t, $,Tj, 0), a 
chosen y-particle, which is close to at time Tj, does not jump on the time interval [r^Tj + e/p]. 

By |16l Prop. 8], the first term in (|88|1 is bounded from above by 1/t 2 for t large enough. For the 
second term we have p(^)|{|t/4J,- ,t}\ < e f° r some C > and all t, while large deviations yield that 
we can find e > such that 

t/2 

fc=l 

for t large enough. From now on we fix such an e. Borel-Cantelli then yields that F^-a.s., G(t) holds 
for all t G N large enough. 

Next observe that by the strong Markov property of X, we can construct a coupling such that on 
the event Efce{L</ 4 J - t-i} ^A(t,x,k,e) > e *> the random variable J Q * £(s, A(s)) ds almost surely dominates 
the sum of i.i.d. random variables (<& e )i<i<et with 

i, e = e/(2p)) =a:= inf _ if = * Vs G [e/(2p), e/p)]) > 0, 



2/^SQr(0) 



= 0) 



I — a: 



qi, e corresponds to the event that given n := inf{j > [t/Al : E 



k=[t/4\ ^A{t,X,k,e) 



i}, X finds a 



y-particle in the £ field which guarantees the event Ait, X, n, e), and then occupies the same position 
as that y-particle on the time interval [r, + e/(2p),Tj + e/p]. Since P^-a.s., G(t) holds for all t G N 
large enough, for such t, we have 



E 



-V 



exp{- 7 j\(s,X(s))d S },~(t) 



< E 



i=l 9i,£ 



( 



-7€/(2 / t>) 



+ 1 - a) 



(89) 



Thus, with dSZD and ([MI) we obtain that P«-a.s., for all i G N large, 



E* [ exp { - 7 J £( s , ds}] < [exp { - 7 ^ ds} , H(t)] + Pff {X G H(t) c ) 



< e 



-S(t+o(t)) 



for some <5 > 0. This establishes the desired result along integer t. Since is monotone in t, we 
deduce that the result holds as stated. I 
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